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We investigate single molecule magnet transistors (SMMTs) with ligands that support transport resonances. 
We find the lowest unoccupied molecular orbitals of Mni2-benzoate SMMs (with and without thiol or methyl- 
sulfide termination) to be on ligands, the highest occupied molecular orbitals being on the Mni2 magnetic core. 
We predict gate controlled switching between Coulomb blockade and coherent resonant tunneling in SMMTs 
based on such SMMs, strong spin filtering by the SMM in both transport regimes, and that if such switching is 
observed then the magnetic easy axis of the SMM is parallel to the direction of the current through the SMM. 
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Charge and spin transport through individual molecules 
bridging a pair of electrodes, and through single-molecule 
transistors that include a third "gate" electrode, have been 
studied intensively for more than a decade.' Because of 
their large magnetic anisotropy barriers and associated stable 
magnetic moments,- single molecule magnets (SMMs) raise 
the possibility of molecular spintronic devices and molecu- 
lar magnetic information storage. Therefore, at present, the 
transport properties of transistors based on individual SMMs 
are attracting considerable experimental^'' and theoretical^'- 
interest.^^ Models based on effective spin-Hamiltonians"*"** and 
density functional theory (DFT)-based calculations'* have 
yielded many important insights into the transport properties 
of these systems. However, in the theoretical transport studies 
to date the organic ligands that surround the magnetic cores of 
the SMMs have received limited attention: The focus has been 
on SMMs in which the lowest unoccupied molecular orbital 
(LUMO) and the highest occupied molecular orbital (HOMO) 
are both on the SMM's magnetic core. Thus the ligands have 
acted as simple tunnel barriers that hinder electron transmis- 
sion between the electrodes and the magnetic core. In the case 
of weak tunneling, the electrostatic energy associated with 
charging of the core states during transport suppresses electri- 
cal conduction at low source-drain and gate voltages.''^ This 
phenomenon is known as Coulomb blockade and has been ob- 
served in SMM transport experiments.^ 

In this paper we explore theoretically the electron and 
spin transport in SMMs in which we predict the ligands to 
play a more active and interesting role: We consider Mni2- 
benzoate [Mni20i2(02CC6H5)i6(H20)4, abbreviated Mni2- 
Ph] SMMs with and without terminating methyl-sulfide or 
thiol groups. For these SMMs we predict the lowest un- 
occupied molecular orbital (LUMO), and orbitals nearby in 
energy, to be on ligands, and the highest occupied molecu- 
lar orbital (HOMO) and orbitals nearby in energy to be on 
the SMM's magnetic core. We also predict that, for certain 
geometries of these SMMs thiol-bonded to gold electrodes, 
some molecular orbitals that are close in energy to the LUMO 
hybridize strongly with the gold electrodes. It follows that 
for these bonding geometries transport via these near-LUMO 
orbitals should not be subject to Coulomb blockade, unlike 
transport via the HOMO. As we explain below SMM transis- 



tors in which such bonding geometries are reaUzed can readily 
be identified experimentally by carrying out measurements of 
the electric current versus source-drain and gate voltage such 
as are routinely performed in single-molecule transistor ex- 
periments today. Furthermore, for SMM transistors identi- 
fied in this way we predict the molecular easy axis to be ap- 
proximately parallel to the direction of current flow. Thus the 
present theory provides a practical way to use standard «o«- 
spin-resolved current measurements to determine the orienta- 
tion of the SMM magnetic easy axis, over which there was no 
control in SMM transistor experiments to date. Moreover, we 
predict these magnetically oriented SMM transistors to be ef- 
fective spin filters at low and moderate source-drain voltages. 

The results that we present are based on the semi-empirical 
extended Hiickel tight-binding model of quantum chemistry 
that we generalize in this paper to include spin polariza- 
tion and spin-orbit coupling. For Mni2 SMMs (both neu- 
tral molecules and negatively charged ions) the present model 
yields results consistent with experiment for the total spin of 
the SMM, for the spins of the individual Mn ions, for the di- 
rection of the magnetic easy axis, for the size of the mag- 
netic anisotropy barrier (MAB), for the size of the molecu- 
lar HOMO-LUMO gap and for the spins of the HOMO and 
LUMO states. '^^ The overall degree of agreement with experi- 
ment obtained with the present model is comparable to or bet- 
ter than that achieved with DFT calculations corrected by in- 
clusion of the adjustable Hubbard U parameter.'' "' However, 
calculations based on the present model are much less com- 
pute intensive than those based on DFT. Thus, we are able to 
study transport in larger molecules than are readily accessible 
to DFT computations and therefore, unlike in previous the- 
oretical studies, to include complete sets of ligands none of 
which have been shortened or replaced by hydrogen atoms. 

Our SMM Hamiltonian is H^^^ = i?™ + iJ^P™ + H^°. 
Here is the extended Hiickel Hamiltonian.'''"* '^ The spin 
Hamiltonian iJ^P'" gives rise to the magnetic polarization of 
the molecule. Spin orbit coupling is described by 

In extended Hiickel theory the basis is a small set of 
Slater-type atomic valence orbitals l^'ia); is the i*"^ 

orbital of the a*'^ atom. The diagonal elements of i?™ 
are the experimentally determined negative valence orbital 
ionization energies e;, (^'ia 1^?™!*^^) = ^f™ia = ^ia- 
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The non-diagonal matrix elements are assumed to be propor- 
tional to the orbital overlaps Dia i'a' = (^ial^i'a')- i-e.. 



and 



EH 



'K , where K is chosen empiri- 

cally for consistency with experimental molecular electronic 
structure data. In our calculations,''' as in Ref. 14, K — 
1.75 + A,L;,,„, - 0.75A4„^^,„, where A,„;,,„, = 

For non-magnetic systems, transport calculations based 
on extended Hiickel theory have yielded elastic'^"''' and 
inelastic"* conductances in agreement with experiment for 
molecules thiol-bonded to gold electrodes and have also 
explained transport phenomena observed in scanning tun- 
neling microscopy (STM) experiments on molecular arrays 
on silicon'^ as well as electroluminescence data"", current- 
voltage characteristics-" and STM images-' of molecules on 
complex substrates. 

The extended Hiickel Hamiltonian i7™ does not describe 
spin polarization whereas in the Mni2 SMMs the four inner 
and eight outer Mn ions are spin polarized with antiparallel 
spins. In our model, H^'p™ addresses this issue. We define 
its matrix elements {isa\H'^'P™\i' s' a') ~ -H^iTaVs'a' between 
valence orbitals i and i' of atoms a and a' with spin s and s' 
by 



{-dinner for inner Mn rf- valence orbitals 
Pouter for outer Mn rf- valence orbitals 
otherwise 



(1) 



Here n is a unit vector aligned with the magnetic moment of 
the SMM, and S is the one-electron spin operator, dinner and 
•Pouter are parameters chosen so that in the Mni2 ground state 
the spin of each inner (outer) Mn is 5'inncr(outor) — ~f (+2)- 
Spin-orbit coupling is also not included in extended Hiickel 
theory. ''''*''^ However, it is responsible for the magnetic 
anisotropy of SMMs. We therefore generalize extended 
Hiickel theory to include spin-orbit coupling by evaluating 
approximately the matrix elements of the spin-orbit coupling 
Hamiltonian H^^ starting from the standard expression-- 
= cr • W{y) X p h/{2mcY where p is the momen- 
tum operator, V{y) is the electron Coulomb potential energy, 
c = (crr, CTy, CT^) and ax,(Jy and az are the Pauli spin ma- 
trices. We approximate ^(r) by a sum of atomic poten- 
tial energies V{r) ~ ^"(i" ~ Tq) where Tq is the posi- 
tion of a^^ atomic nucleus. Noting that the spin-orbit cou- 
pling arises mainly from the atomic cores where the poten- 
tial energy Va{r — r^) is approximately spherically symmet- 
ric, yields H^^ ^ ^_^^^V|iJ^S . ^here 

S = her /2, Lq, = (r — r„) X p is the orbital angular momen- 
tum operator with respect to the position of the a*'^ nucleus 
and the sum is over all atoms a. Evaluating the matrix ele- 
ments of between valence orbitals i and i' of atoms a 
and a' with spin s and s' we then find 



so 



(2) 



The first term on the right-hand side is the intra-atomic contri- 
bution, the remaining terms are the inter-atomic contribution 



^so 



X {R^ 



dV(\y-r^\) (3) 

-^\Rc.,u)oi,i., 



A d(\ 



where the atomic orbital wave function "^ia has been ex- 
pressed as the product of a radial wave function Ra,u 
and directed atomic orbital \a,li,di,s). Here li is 
the angular momentum quantum number and di may be 



^ ; Px 7 Py 1 Vz 7 dxy 7 dx 



depending on the value of li. In Eq. 



(3) the matrix elements of S ■ are evaluated as in Ref. 23 
while the radial integrals are the spin-orbit coupling constants. 

We estimate the magnetic anisotropy energy of the SMM 
from the total ground state energy expression i^totai = 
Y.iEi-Y.i{4>i 0i) /2 where Ei and are eigenen- 

ergies and eigenstates of i/^MM ^j^^ j-j^g summations are over 
all occupied states. Since i7*^P™ represents electron-electron 
interactions at a mean field level, the second summation on 
the right hand side is required to avoid double counting the 
corresponding interaction energy. 

Experimental estimates-"* of the spin-orbit coupling con- 
stant for Mn (i-orbitals have been in the range 0.023-0.05 1 e V, 
while theoretical estimates-'' have been in the range 0.038- 
0.055 eV. In this paper for Mn atoms we use the value 
0.036 eV, which is consistent with the experimental and theo- 
retical values. We find that the spin-orbit coupling constants 
of the other atoms do not affect the SMMs' properties sig- 
nificantly and that the intra- and inter-atomic terms in Eq. 
(2) make contributions of the same order of magnitude to the 
magnetic anisotropy barriers of Mni2 SMMs. In this work the 
values of Hf^^,^, and Dia^ua' that enter the extended Hiickel 
model, Eqs. (1) and (2), were adopted without change from 
Refs. 14 and 15. The molecular geometries that we stud- 
ied were based on the experimentally measured geometry^^ 
of Mni2-Ph modified as necessary by adding thiol or methyl 
sulfide groups to the ligands. Thus the only free parameters 
in the present theory are the Aia of Eq. (1) that control the 
spin polarizations of the Mn atoms; we chose dinner = 3.0 
eV and Pouter = —3.5 eV. For these parameter values in the 
Mni2-Ph ground state we find the inner and outer Mn ions to 
have spins -1.6 and 1.99, the SMM to have a total spin of 10 
and the calculated MAB to be 2.50 meV, all consistent with 
experiment. -^^'-^ 

The calculated densities of states projected on the inner and 
outer Mn atoms, and carbon atoms of Mni2-Ph are shown in 
Fig. 1. The HOMO and nearby levels are on the outer Mn 
and are filled with spin-up electrons (parallel to the total spin). 



consistent with S„ 



~2. The occupied inner Mn states are 



filled with spin-down electrons consistent with 5innei = — f ■ 
The carbon atoms are weakly spin polarized. The calculated 
HOMO-LUMO energy gap is '^0.7eV. (Note that experimen- 
tal and theoretical estimates for the Mni2 family range from a 
few tenths of an electron Volt to more than 1 eV. '"''- -**) As is 
seen in Fig. 1, our calculation predicts the LUMO of Mni2- 
Ph to be mainly on the carbon atoms of the ligands. There 
is at this time no direct experimental evidence as to whether 
the LUMO of Mni2-Ph is located on the ligands as in Fig. 1 
or on the Mni2 core of the molecule. However, experimental 
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Figure 1. Projected density of states for majority (spin-up) and mi- 
nority (spin-down) electrons on Mn and carbon atoms for tiie isolated 
Mni2-Ph molecule without thiol or methylsulfide end groups. 



measurements of the MAB carried out on the neutral-^ and 
negatively charged"'^ Mni2-Ph species yielded values of 3.3 
and 2.41 meV, respectively, a difference of only 27%. Since 
the main source of the MAB of Mn 12 is the Jahn-Teller distor- 
tion of the Mni2 core of the molecule^ this insensitivity of the 
MAB to the oxidation state of the molecule is consistent with 
the added electron of the negatively charged species residing 
mainly on the ligands rather than on the Mni2 core, as one 
might expect if the LUMO of the neutral molecule is on the 
ligands as in Fig. 1 . Our calculations also predict only a small 
change in the MAB when an electron is added to the neutral 
Mni2-Ph molecule, since in our model the added electron lo- 
cates primarily on the ligands rather than the magnetic molec- 
ular core. Our prediction that the LUMO is on the ligands(i.e. 
the benzoate groups, including both their carbon and oxygen 
atoms) and not on the Mni2 core is also consistent with the 
large 3.3-3.6 eV electron affinity of the benzoate species. 

The ligands of SMMs that are studied in transport experi- 
ments with gold electrodes are normally thiolated and we shall 
therefore focus our attention on Mni2-Ph-Th, i.e., Mni2-Ph 
terminated with methyl sulfide (SCH3) groups, the methyl be- 



ing displaced by gold when the molecule bonds to the con- 
tacts. The calculated densities of states of Mni2-Ph-Th (and 
of Mni2-Ph with thiol end groups) projected on the Mn and C 
atoms are similar to those shown in Fig. 1 for Mni2-Ph. The 
locations and spins of the HOMO and LUMO and of orbitals 
close to these in energy are also similar for Mni2-Ph with and 
without thiol or methyl sulfide end groups. The calculated 
HOMO, LUMO and another representative molecular orbital 
close in energy to the LUMO are shown in Fig. 2 for Mni2- 
Ph-Th. The HOMO is located on the magnetic core of the 
molecule. However, the LUMO and molecular orbitals close 
in energy to the LUMO are located on ligands, and specifi- 
cally on those ligands that are oriented approximately paral- 
lel to the magnetic easy axis of the molecule which points in 
the z-direction in Fig. 2. This finding has important impli- 
cations for electron and spin transport in Mni2-Ph-Th-based 
SMM transistors: The HOMO and molecular orbitals close in 
energy to the HOMO have very little overlap with any of the 
ligands. Therefore all of the ligands that couple the molecule 
to the contacts act as strong tunnel barriers for transport me- 
diated by the HOMO and molecular orbitals nearby in energy. 
For this reason transport via the HOMO and molecular or- 
bitals close in energy to the HOMO is predicted to display 
the classic signature of Coulomb blockade. By contrast if any 
gold contact bonds to a ligand on which certain molecular or- 
bitals close in energy to the LUMO (for example, that in Fig. 
2(c)) have a strong presence, those molecular orbitals will hy- 
bridize strongly with the gold contact and therefore transport 
via that molecular orbital (or orbitals) will not be subject to 
Coulomb blockade. Furthermore, because these orbitals oc- 
cupy ligands that are oriented approximately parallel to the 
molecular easy axis as in Fig. 2, if such transport that is not 
subject to Coulomb blockade is observed experimentally, then 
the molecule must be oriented relative to the gold contacts 
in such a way that the magnetic easy axis is approximately 
parallel to the direction of cuiTent flow through the molecule. 
We note that although there have been theoretical suggestions 
previously of possible ways to determine the orientation of 
the easy axis in a SMM transistor,**''- these suggestions have 
been difficult to implement in practice and no experimental 
control over the orientation of the easy axis in SMM transi- 
tors has been achieved experimentally to date. The present 
theory is much more promising in this regard since the exper- 
imental observation of the presence or absence of Coulomb 




Figure 2. Wave functions for Mni2-Ph-Th (a) spin-up HOMO (b) 
spin-down LUMO (c) spin-up orbital near LUMO in energy. The 
magnetic easy axis and total molecular spin are parallel to the z- 
axis. Although the LUMO and near LUMO orbitals are mainly on 
the ligands they are spin polarized due to their small but non-zero 
overlaps with the magnetic core of the molecule that are visible in 
(b) and (c). 




0.2 0.4 0.6 0.2 04 0.6 
Bias voltage (eV) Bias voltage (eV) 



Figure 3. Calculated spin resolved current parallel (||) and perpen- 
dicular (±) to the easy axis for (a) positive and (b) negative gate 
voltage vs. bias voltage at zero temperature. The gate potential at 
the molecule = -1-0.2 and -0.2V in (a) and (b), respectively. 
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blockade in single-molecule transistors is currently carried 
out routinely, and we predict that if Coulomb blockade is ob- 
served at negative gate voltages (transport via the HOMO and 
nearby states) but not at positive gate voltages (transport via 
the LUMO or nearby states) then the easy axis is approxi- 
mately parallel to the direction of current flow. 

Our transport calculations based on Landauer theory and 
the Lippmann-Schwinger equation"*- predict Mni2-Ph-Th- 
based SMM transistors with gold contacts to be effective spin 
filters at low source-drain bias for both positive and nega- 
tive gate voltages. Representative results for a Mni2-Ph-Th 
molecule bonded to gold contacts via ligands on which near- 
LUMO orbitals have a strong presence (and therefore the cur- 
rent through the SMM is roughly parallel to the magnetic easy 
axis) are shown in Fig. 3. Here mainly spin-up electrons are 
transmitted through the SMM at low bias for both positive 
(Fig. 3(a)) and negative (Fig. 3(b)) gate voltages. The gradual 
rise of the current with bias voltage (from ^0.2 to ^0.67V) 
in Fig. 3(a) is a direct manifestation of the large broadening of 
the near-LUMO molecular orbitals responsible for transport 
that is due to the strong hybridization of those orbitals with 
the gold contacts that also suppresses Coulomb blockade for 
positive gate bias. By contrast, the abrupt step-like rise to 
much lower values of the current for negative gate voltages 
(Fig. 3(b)) is due to the near-HOMO molecular orbitals being 
very weakly coupled to the gold contacts and therefore be- 
ing only very weakly broadened, they are subject to Coulomb 
blockade. For molecules bonded to gold electrodes via lig- 
ands that are roughly perpendicular to the easy axis, neither 
near-LUMO nor near-HOMO orbitals have significant pres- 
ence on the ligands, see Fig. 2. Thus, for geometries with a 
current through the SMM that is roughly perpendicular to the 
magnetic easy axis, the ligands act as strong tunnel barriers 
for both positive and negative gate voltages (Fig. 3). Coulomb 
blockade is therefore predicted for both signs of the gate volt- 
age. 

In conclusion, previous studies of transport in single 
molecule magnets have considered the situation where the 
HOMO and LUMO both reside on the magnetic core of the 
molecule. Here we have proposed that this need not al- 
ways be the case, Mni2 benzoate with and without termi- 



nating methyl-sulfide or thiol groups being a possible exam- 
ple. We have predicted that for these systems the LUMO 
and molecular orbitals close in energy to the LUMO re- 
side on the organic ligands that are oriented approximately 
parallel to the magnetic easy axis of the molecule and that 
when the molecule bonds via these ligands to gold electrodes 
in a single-molecule transistor, transport via some of the 
near-LUMO orbitals should not be in the Coulomb blockade 
regime. For other orientations of the molecule transport via 
the LUMO and near-LUMO orbitals is predicted be in the 
Coulomb blockade regime, as is transport via the HOMO and 
near-HOMO orbitals for all molecular orientations. This ef- 
fect should make it possible to study experimentally the trans- 
port in single molecule magnet transistors that behave as spin 
filters and in which the orientation of the magnetic easy axis 
relative to the electrodes is known. In single-molecule magnet 
transistor experiments to date the orientation of the magnetic 
easy axis has not been determined although it controls the spin 
polarization of the current in such devices. 

While this work has focused on Mni2 benzoate and its 
derivatives, we expect other single molecule magnets with 
LUMO and/or HOMO states located on the ligands to exist as 
well due to the small HOMO-LUMO gaps exhibited by a va- 
riety of organic molecules that may be chosen as ligands. For 
example it is well established that polyacetylene and polythio- 
phene have HOMO-LUMO gaps of L4 eV" and 0.85 eV^^^ 
respectively. Both of these are smaller than the experimentally 
measured energy gap (^ L8 eV)^^'* between the highest oc- 
cupied and lowest unoccupied orbitals localized on the cores 
of the Mni2 molecules. It is reasonable to expect oligomers 
(consisting of several monomers) derived from these poly- 
mers to have similarly small HOMO-LUMO gaps. If such 
oligomers are used as ligands for Mni2 SMMs contacted with 
gold electrodes, the gold Fermi level is expected to lie within 
the oligomer HOMO-LUMO gap and also within the HOMO- 
LUMO gap of the Mni2 core. Therefore either the HOMO or 
the LUMO (or both) of the Mni2 SMM with such ligands is 
expected to lie on the ligands. 
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